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Abstract. This talk will focus on recent results by the MILC collaboration from simulations 
of light hadrons in 2+1 flavor lattice QCD. We have achieved high precision results in the 
pseudoscalar sector, including masses and decay constants, plus quark masses and Gasser- 
Leutwyler parameters from well controlled chiral perturbation theory fits to our data. We 
also show spectroscopy results for vector mesons and baryons. 



1. Introduction 

The MILC Collaboration has been engaged in the study of QCD with three flavors of improved 
staggered sea quarks for some time. Our goals are to determine the hadron mass spectrum, 
the properties of light pseudoscalar mesons, the topological susceptibility, the decay properties 
of B and D mesons, and the behavior of strongly interacting matter at high temperatures, all 
in as realistic simulations as possible. In this talk I will concentrate on the light pseudoscalar 
sector and the light hadron spectroscopy [2J |3] . 

While we have made simulations with lighter u and d quark masses (with m u = = mi) 
than have been reached before, to simulate at the physical mass value has not been possible 
so far even with the improved "asqtad" quarks pQ. We have created gauge field ensembles at 
two lattice spacings, a ~ 0.125 and ~ 0.09 fm, with the strange quark mass m s fixed near its 
physical value and several light quark masses, as low as mi = 0.1 m s . This allows controlled chiral 
extrapolations to the physical light quark mass, and from the two lattice spacings, extrapolations 
to the continuum limit, a — ► 0. The physical spatial size of our gauge field ensembles is L > 2.5 
fm. At one set of simulation parameters we created ensembles with two physical sizes, L = 2.5 
and 3.5 fm, to check for finite volume effects. These were found to be small. 

All gauge field ensembles created by the MILC collaboration are made available to other 
researchers at the "NERSC Gauge Connection". They have been used, for example, by the 
Fermilab, HPQCD and UKQCD collaborations for the study of heavy quarkonia and heavy- 
light mesons. A comparison of quantities that can be computed on the lattice with small 
systematic uncertainties and that are well known experimentally shows agreement within errors 

of 1-3% m. 



2. The light pseudoscalar sector 

In the pseudoscalar sector we have the most accurate data with many partially quenched 
measurements, i.e., measurements of pseudoscalar masses and decay constants with valence 
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quark masses different from the sea quark masses. Fits, GVGI1 clt cl single lattice spacing, to 
continuum partially quenched chiral perturbation theory (XPT) do not work Lattice effects, 
in particular the 0{a 2 s a 2 ) taste breaking effects of the improved staggered fermions, need to be 
taken into account. 

The staggered XPT (SXPT) formalism for our case of 2+1 sea quark flavors with partially 
quenched measurements has been worked out in detail in [5]. All our data for f w and m 2 are fit 
to the SXPT form simultaneously, taking account of the taste violations. The continuum XPT 
parameters were allowed to have (D(a s a 2 ) terms. After the fit, all the terms taking account of 
the finite lattice spacing effects were set to zero to take the continuum limit. Small finite volume 
errors (< 1.5%) were corrected within the XPT framework. Details can be found in [2]. Sample 
plots for "full QCD" , i. e. with m va i = m sea are shown in Fig. ^ 
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Figure 1. Illustration of SXPT fit showing full QCD points for f w (left) and m 2 /(m x + m y ) 
(right). Also shown are the fits extrapolated to the continuum limit, and after correcting from 
the simulation strange quark mass to its physical value. 



The strange quark mass in the simulations, our best guess beforehand and denoted by m' s in 
Fig. E turned out a little larger than the physical strange quark mass, m s . We can correct for 
this within the chiral fits. The corresponding curves are also shown in Fig. ^ 

To find the quark masses, we must extrapolate to the physical meson masses. Electromagnetic 
and isospin- violating effects {m u ^ m^) are important. Denoting by and m^- the masses 
with EM effects turned off and m u = = m, i.e. what is done in our simulations, we have 
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where the superscript QCD indicates the masses with EM effects turned off. A# = is "Dashen's 
theorem." Continuum considerations suggest that A^; ~ 1. 



With the considerations in eq. Q we found, in collaboration with the HPQCD and UKQCD 
groups, using a one-loop mass renormalization constant, [H] 

mf (2 GeV) = 76(0) (3) (7) (0) MeV , 

m m (2 GeV) = 2.8(0)(1)(3)(0) MeV , (2) 
m s /rh = 27.4(1) (4) (0)(1) . 

Here, the errors are from statistics, simulation, perturbation theory, and EM effects (obtained 
from varying between and 2), respectively. 
For the pseudoscalar decay constants we obtain 

f n = 129.5 ± 0.9 ± 3.5 MeV , 

fx = 156.6 ± 1.0 ± 3.6 MeV , (3) 
fx/U = 1.210(4)(13) , 

where the first error is statistical and the second is systematic. More details, including our 
results for the Gasser-Leutwyler chiral parameters and the ratio m„/m,j can be found in 

Marciano has pointed out [7j that an accurate determination of fx / fn can be used to calculate 
the CKM matrix element V us . With the above result we obtain V us = 0.2219(26). The error, 
already comparable to that from the standard method using decays, is dominated by the 
lattice error on fx/ fn- With the coming simulations, V us will be known more precisely. 



3. Light hadron spectroscopy 
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Figure 2. Masses of the groundstate and first excited pseudoscalar mesons for pion like states 
(left) and kaon like states (right). 



In the pseudoscalar sector the statistical errors on the correlation functions were small enough 
that we could get some signal for excited states, as shown in Fig. [2J It is noteworthy that in 
kaon correlation functions we could observe an opposite parity state - correlation functions with 
staggered fermions generically contain opposite parity states - with energy close to m n + mx, 
i.e. the expected decay products of the opposite parity excited meson. No such state was found 
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Figure 3. Vector meson masses. 
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Figure 4. Nucleon masses. 
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in matched quenched simulations - absence of virtual quark loops in quenched QCD prevents 
emergence of the ll ir + K" intermediate state. 

Our vector meson and nucleon masses are shown in Figs. 01 and In the latter we also 
show some possible chiral extrapolations, after taking the continuum limit at fixed mi/m s . For 
further details and more spectroscopy results see [2]. 

We end by showing a summary of our current spectroscopy results in Fig. El 
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Figure 5. Comparison of our hadron spec- 
troscopy, from a crude chiral and contin- 
uum extrapolation, with experiment, tt and 
K were used to set the light and strange 
quark masses. The upsilon and charmonium 
columns are differences from the ground state 
(IS) masses, from work of the HPQCD and 
Fermilab groups [UIB1- The T IP-IS splitting 
was used to fix the lattice spacing. 
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